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Even in a topologically perfect crystal, a moving twin wall will experience forces due to the discrete nature 
of the lattice. The potential energy landscape can be described in terms of one of two parameters: the 
Peierls energy, which is the activation energy for domain wall motion in a perfect crystal; and the Peierls 
stress, the maximum pinning stress that the potential can exert. We investigate these parameters in a one 
order parameter discrete Landau-Ginzburg model and a classical potential model of the ferroclastic perovskite 
CaTiOs- Using the one order parameter model we show that the Peierls energy scales with the barrier height 
of the Landau double well potential and calculate its dependence on the width of the wall numerically. In 
CaTiOs we calculate the Peierls energy and stress indirectly from the one order parameter model and directly 
from the interatomic force field. Despite the simplicity of the one order parameter model, its predictions of 
the activation energy are in good agreement with calculated values. 

FAGS numbers: 62.20.Dc, 61.72.Mm 
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I. INTRODUCTION 

The motion of ferroclastic or ferroelectric-ferroelastic 
twin walls plays a significant role in determining the elas- 
tic, dielectric, piezoelectric and ferroelectric properties 
of a number of materials of scientific and technological 
interest The lower mantle of the Earth is known to 
consist mainly of magnesium silicate perovskite, a ferroc- 
lastic polymorph of MgSiOa. Recent work has explored 
the possibility that the seismic properties of the lower 
mantle, such as attenuation, can be explained in terms 
of the elastic response of domain walls. The large piezo- 
electric and dielectric coefficients of barium titanate and 
lead titanate have been shown to have significant con- 
tributions from the motion of twin walls. Finally, ferro- 
electric switching, which is currently being exploited for 
computer random access memory applications, is known 
to be entirely due to the motion of twin walls. 

To understand the properties of these materials, and 
the systems in which they are found, we must under- 
stand the factors which affect the motion of twin walls. 
Unlike magnetic domain walls, with widths of lOO's of 
nanometres, ferroclastic and ferroelectric walls are atom- 
istically thin, with wall widths of the order of the unit cell 
parameterii We must understand their behaviour from 
an atomic perspective. In essence this requires an un- 
derstanding of the energy landscape through which twin 
walls move. 

A schematic of this energy landscape is shown in Fig.[T] 
As a wall moves through the crystal it experiences a po- 
tential which oscillates between i^minimum at local minima 
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Figure 1. The energy landscape experienced by a domain 
wall moving through a crystal containing a point defect. The 
bold line shows the potential experienced by the wall while 
the thin lines show the local equilibrium energy of the wall 
iSminimum and the Saddle point energy of the transition state 
from one local minimum to another -Bsaddio- Far from a defect 
the difference between these two energies, i.e. the activation 
energy needed for the wall to move is called the Peierls energy 

-Epcicrls- 



and i?saddic at transition states (saddle points) between 
two minima. If the motion of the wall is a thermally ac- 
tivated process then the mobility of the wall depends on 
the difference between these two energies. 
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As shown in the figure there are two contributions to 
the potential energy landscape. These can be labelled 
intrinsic and extrinsic. The intrinsic contribution to the 
energy landscape is present even in a chemically and 
topologically perfect crystal and is due to the periodic- 
ity of the lattice. This contribution is parameterised by 
the Peierls energy -Epdcris , which is the activation energy 
for twin wall motion far from point defects. The second, 
extrinsic contribution to the energy landscape is due to 



defects in the perfect lattice such as vacancies, impurity 
atoms, dislocations and other twin walls. 

There have been attempts to understand this energy 
landscape using both experimental and simulation meth- 
ods. Experimentally, it is clear that extrinsic pinning 
due to point defects is far more significant than lattice 
pinning,— i^i^ which is often too small to detect, except 
by very sensitive metliodsi^ Simulations of oxygen va- 
cancies in ferroelastic calcium titanateS and ferroelectric 
lead zirconatai^ have already been carried out, show- 
ing that an oxygen vacancy has an energy approximately 
1 eV lower in the wall than in the bulk. These simula- 
tions provide information about -Eminimum, but in order 
to complete the picture, calculations of -Esaddic arc also 
necessary. 

In this work we only consider intrinsic pinning. We 
investigate pinning in a one order parameter discrete 
Landau-Ginzburg model, building on previous work by 
Ishibashiii and Combs and Yipi^ii^. Then we investi- 
gate intrinsic pinning in an empirical potential model 
of orthorhombic calcium titanatc CaTiOs developed by 
Calleja et ali^ We compare the predictions of the one or- 
der parameter model with the results of a transition state 
calculation and show a good agreement between the two 
values. 



II. A ONE PARAMETER MODEL 

In this section wc investigate the intrinsic pinning of 
domain walls in the discrete Landau-Ginzburg or 0^ 
models^ Intrinsic pinning in this model has been inves- 
tigated by Ishibashi and Combs and Yipi^i"— We re- 
port the results of a numerical calculation of the pinning 
energy showing, in agreement with previous work, that 
when the wall width is two lattice spacings the activation 
energy is practically zero. 

The most successful theoretical tool for describing 
phase transitions in ferroelectric and ferroelastic materi- 
als is Landau-Ginzburg theory J^ii^ Through the Landau- 
Ginzburg free energy the theory provides a framework 
which can be used to predict both macroscopic be- 
haviour, such as the specific heat capacityi^ and elastic 
constantsii of a material going through a phase transi- 
tion, and microstructural details, such as the structure 
of domain walls J^ii^ 

Usually a continuum formulation of the Landau- 
Ginzburg free energy is used, in which the discrete nature 
of the lattice is neglected i^ii^ This approach has been 
very successful even in predicting the structure of twin 
walls, where the continuum approximation might be ex- 
pected to break down. It is not possible to calculate the 
Peierls energy within the continuum limit and so we use 
a discrete form of the Landau-Ginzburg energy. 



The first term is a double well potential, where AE is 
height of the barrier between the two walls. The second 
term is the discrete analogue of the Ginzburg term, w 
is the wall width and a is the lattice spacing. These 
quantities are illustrated in Fig. [2] 
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Figure 2. Dimensional quantities of the general model, 
(a) The Landau double well potential is characterised by an 
energy barrier AE. (b) The discrete nature of the lattice 
means the model contains two lengthscales: a the lattice pa- 
rameter and w the wall width. 

Dimensional analysis tells us that the Peierls energy 
-Bpoioris must be given by 



£^Poioris = AEf (w/a) 



(3) 



where / is to be determined. It is easy to deduce the 
limiting values of f{x) in the cases when x is very small 
or very large. 

In the case u; = the free energy of the discrete 
Landau-Ginzburg model is 



f = J2ae{q^-i) 



(4) 



The system consists of a collection of independent order 
parameters Qi moving in double well potentials. The do- 
main wall moves when one order parameter flips from one 
state to another. The activation energy for this process 
is AE and thus 



lim/(a;) = l 



(5) 



If w is very large then the discrete nature of the lattice is 
irrelevant and a continuum approximation may be used. 
In the continuum theory the energy of a domain wall is 
independent of its position, and thus there is no activa- 
tion energy for domain wall motion 



lim / I 
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For intermediate cases we calculate f{x) numerically. 
We calculate the energy of a single wall in a 200 site 
system, both without constraints (i?minimum) and with 
the constraint that Qioo = (£'saddio)- The difference 
between the two energies divided by AE gives us f{x), 
shown in Fig [3l The figure shows that when the wall 
width is twice the lattice parameter the Peierls energy is 
already practically zero. 

In a crystal the width of a twin wall can be affected 
by two parameters: temperature and velocity. Elemen- 
tary Landau-Ginzburg theory predicts that the twin wall 
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Figure 3. Dependence of the activation energy on the ratio 
between wall width and lattice parameter. 



width should diverge as T approaches Tc.— This predic- 
tion has been confirmed experimentally^^ As the temper- 
ature approaches the activation energy for wall motion 
will decrease. The second factor which can affect the wall 
width is the velocity of the wall. If the speed of the wall 
V approaches the velocity of sound in the material c then 
the width of the wall is 'Lorentz contracted' by a factor 
of 



(7) 



As a wall accelerates the forces it experiences due to the 
lattice potential increase. 



III. TWIN WALL MOTION IN CATIO3 

To test the validity of the above approach we com- 
pared the value of the activation energy calculated by the 
method described above with a direct transition state en- 
ergy calculation. We used an empirical potentials model 
to investigate twin walls in CaTiOa. We calculate the 
structure of the twin walls of the system and calculate 
the Peierls energy and stress. 

CaTiOa is a ferroelastic, but not ferroelectric, per- 
ovskite. The crystal structure consists of corner linked 
TiOg octahedra with Ca atoms distributed between the 
octahedra. At high temperatures the crystal structure is 
cubic but at room temperature the crystal structure is 
orthorhombic, with a space group of Pbnm and a Glazer 
octahedral tilt system of a~a~c~^^^^ The crystal struc- 
ture and the coordinate system used in this work is shown 
in Fig. m When measurements of wall widths are given 
below they are given in units of the pseudocubic unit cell, 
containing a single formula unit. 



A. Structure of static twin walls. 

In this work we consider a ferroelastic wall perpen- 
dicular to the X-axis. The structure of the wall can 
be described in terms of order parameters and strains. 
The Glazer tilt system allows us to define order parame- 
ters {Qx^QyiQ z) associated with rotations of octahedra 
about the x-, y-, and z-axes {0x,&y,&z)- If the position 



Figure 4. Crystal structure of CaTiOa, showing the coordi- 
nate system used in this work. Ca atoms are shown in black, 
the octahedra have O atoms at their vertices and Ti atoms at 
the centres. 



of an octahedron in the crystal is labelled by integers 
(ixi iyjiz), then the order parameters are defined by 



Q,(_l)»==+»« 



(8) 
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The compatibility conditions limit the strains which can 
vary across an interface. For an interface perpendicular 
to the a;-axis only the strains exx , ^xy and ezx can be non- 
zero. Furthermore the symmetry of the crystal constrains 
the strain £zx to be zero. The strain £xy is the ferroelastic 
strain. This changes sign across a ferroelastic wall. The 
strain Cxx is a secondary strain, which only takes non-zero 
values within the wall. 

Calleja et ali^ developed an empirical potential set 
for this mineral and investigated the interaction between 
oxygen vacancies and twin walls in a configuration con- 
taining 26 X 10 X 6 cells. The authors simulated a single 
domain structure and then rotated part of their config- 
uration through 90° to generate twin walls. This proce- 
dure generates an interface consisting of the combination 
of a ferroelastic twin wall (with an order parameter Qy) 
with an antiphase boundary (with an order parameter 
Qz)- These two types of wall can exist independently so 
in this work we consider simple ferroelastic twin walls in a 
system of 14 X 6 X 4 octahedra implemented in DL_POLY— 
using Calleja et al.'s potential set. Periodic boundary 
conditions make it impossible to simulate a single twin 
wall so instead we simulate a system with two walls. The 
order parameters and strains in the walls relaxed at abso- 
lute zero are shown in Figs. [5] and [6l The order parameter 
Qy and the shear strain exy change sign across the walls. 
Fitting Qy to a hyperbolic tangent profile gives a wall 
width w = 1.3a. Qx, Qz, and exx show anomalies across 
the wall. This is the behaviour expected from secondary 
order parameters 



B. Activation energy for twin wall motion in CaTiOs. 

We compare two methods of calculating the Peierls en- 
ergy Epcicris- The first method is an indirect calculation 
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Figure 5. The behaviour of the order parameters across the 
wall. The order parameters are described in equation [S] Qy 
shows a hyperbolic tangent variation across the wall, and 
and Qy show anomalies at the wall. 



Indirect Calculation 
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0.116 Jm~^ 


AE 


0.034 J 


w/a 


1.13 


f{w/a) 


0.016 


-^Peierls 


0.530 mJm"^ 


^ Poicrls 


4.35 MPa 


Direct Calculation 


-£'Peierls 


0.313 mJm~^ 


^ Peierls 


2.57 MPa 



Table I. Indirect and direct calculation of the activation en- 
ergy and Peierls stress. The results of the calculations are in 
good agreement. 
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Figure 6. Strain behaviour across the twin wall. The fer- 
roelastic strain e^y follows the order parameter Qy, while the 
secondary strain exx shows an anomaly within the wall. 



sary because DL_POLY cannot directly resolve the energy 
differences involved.) Again the results of these calcula- 
tions are summarised in Tab. HI 

The agreement betvifeen our two results is very good — 
less than a factor of two — especially given the small value 
of the Peierls energy compared with the interaction en- 
ergy of a twin wall with an oxygen vacancy, which, as 
noted above, is of the order of 1 eV. The residual dis- 
crepancy may be due to the complexity of the system. 
Equation [3] was developed for a domain wall which can 
be described by a single order parameter. As shown in 
Figs. [5] and [6] at least five parameters show anomalies 
across the wall. The energies of these anomalies may 
lead to an overestimation of AE calculated from the in- 
terfacial energy 7 of the wall. 



using equation [3l Secondly we perform a direct calcula- 
tion of the transition state energy. These two energies are 
in good agreement. We also calculate the Peierls stress 
fpoioris which is the maximum restoring stress the wall 
can exert. 

To calculate the Peierls energy iJpoioris using equa- 
tion[3l we first calculate AE from the interfacial energy 7 
of the twin walls. From AE and the value of the function 
/ for w/a = 1.3 we can calculate iJpoioris- Finally, assum- 
ing a sinusoidal variation of the energy as the wall moves 
through the lattice we can calculate the Peierls stress 
cpeieris- The stagcs of the calculation are summarised in 
Tab. II] 

To calculate the Peierls energy directly we perform a 
transition state calculation on our system. We start with 
initial and final states which differ in that the two twin 
walls of the system are each translated by one lattice pa- 
rameter in the same direction. The difference in atomic 
coordinates between these two states defines the reaction 
coordinate. We moved the system from the initial to the 
final state in 100 steps, relaxing all the degrees of freedom 
perpendicular to the reaction coordinate. At each step 
we calculated the force along the reaction coordinate and, 
by numerical integration of the work done by that force, 
the Peierls energy and stress. (This approach was neces- 



C. Simulation of a moving domain wall 

Our results suggest that if a pressure greater than 
cpcicris ^ 3 MPa is applied to a twin wall it will move 
freely, rather than as a thermally activated process. In 
this section we demonstrate that this is the case by molec- 
ular dynamics simulation. Working in an NVT ensemble 
we shear the system to generate a force on the walls and 
observe their motion. 

In order to calculate the force on the wall generated 
by a shear stress we need to calculate the Eshelby force 
on the walli^ The stress on a wall crwaii generated by an 
externally applied shear stress a^y is given by 

CWall = 'iCTxyCxy (11) 

where e^y is the spontaneous strain of the transition. For 
CaTiOs €xy^4:X 10"^ (see Fig.^. 

We started with a configuration of 26 x 10 x 6 octahedra 
from the simulation of Calleja et al, containing, as noted 
above, both ferroelastic walls and antiphase boundaries. 
On annealing at 10 K using DL_POLY the antiphase bound- 
aries spontaneously moved together and annihilated each 
other, leaving only the ferroelastic walls. DL_POLY does 
not allow the direct imposition of a constant shear stress 
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Timestep 1.0 fs 

Thermostat relaxation time 0.5 ps 
Simulation duration 10.0 ps 

Initial shear stress on crystal 6.0 GPa 
Initial pressure on walls 100.0 MPa 

Peierls stress 3.0 MPa 



Table II. Parameters used in the simulation of a moving twin 
wall. The stresses acting on the system and the wall are only 
initial values. As the walls move in response to the forces 
these stresses will relax. 
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Figure 7. Observed motion of the two twin walls of the 
simulated system. 



SO instead we sheared the whole system (both coordi- 
nates and velocities) through an angle of 0.3°, generating 
an initial shear stress. (The NVT ensemble prevents the 
relaxation of this stress by a macroscopic shear of the sys- 
tem.) The Eshelby force on the wall exceeded the Peierls 
stress and so motion of the wall was observed. The pa- 
rameters of the simulation are summarised in Tab. [ill 
In response to these forces the walls move as shown in 
Fig. [71 Initially the walls accelerate because the pressure 
acting on them is higher than the Peierls stress. The 
walls traverse several unit cells and then decelerate, as 
the stress acting on them decreases. 



IV. CONCLUSIONS 

A complete picture of twin wall motion in ferroelastic 
and ferroelectric materials would shed light on questions 



such as the fatigue problem in ferroelectric memories 
and the contribution of twin wall motion to the seismic 
properties of the Earth's lower mantle. Such a picture 
requires an understanding of the energy landscape 
through which the twin wall moves in the presence and 
absence of point defects. We have shown that the Peierls 
energy and stress of a fcrroic material can be accurately 
estimated using an indirect approach by mapping the 
system on to a one order parameter model. 
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